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The sixth-order electron-loop vacuum-polarization contri- 
bution to the 2Px/2 — 2Si/2 Lamb shift of the muonic hydro- 
gen (n~p + bound state) is evaluated numerically. Our result 
is 0.007608 (1) meV. This eliminates the largest theoretical 
uncertainty. Combined with the proposed precision measure- 
ment of the Lamb shift it will lead to a precise determination 
of the proton charge radius. 

PACS numbers: 36.10.Dr, 12.20.Ds, 31.30.Jv, 06.20.Jr 



The muonic hydrogen, the p~p + bound state, differs 
from the ordinary hydrogen atom in two important re- 
spects. One is that the vacuum-polarization effect is 
much more important than other radiative corrections. 
The other is that it is more sensitive to the hadronic 
structure of the proton. Thus it provides a means of 
testing aspects of QED significantly different from those 
of the hydrogen atom. 

The muonic hydrogen has a long-lived 25 meta-stablc 
state. This makes it possible to measure the 2P 1 / 2 — 2S , 1 / 2 
Lamb shift to about 10 ppm level using the phase- 
space compressed muon beam technique |lj]. At present, 
however, theoretical precision is limited to about 50 
ppm. This uncertainty comes mainly from the unknown 
contribution AE^ of the sixth-order electron vacuum- 
polarization effect H . 

In this paper we report the result of our evaluation of 
A^ 6 ) . Our result is 



AE (6) = 0.120 045 (12) mJZa) 2 (- 
= 0.007 608 (1) meV, 



(1] 



where Z = 1 for the proton and m r is the reduced mass 
of the p~p + system: j| 



m,,m 



rn. 



v 



= 94.964 485 (28) MeV, 

up T Tflp 

m M = 105.658 389(34) MeV, 



in. 



rn h 



938.272 31(28) MeV 



(2) 
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We have also evaluated the main part of AEW using 
the Pade approximation of vacuum-polarization function 
H . The result ( p4[ ) is in good agreement with the direct 
calculation (|23]). 

The contribution to the 2Pi/ 2 — 2S\/ 2 Lamb Shift of 
the muonic hydrogen due to the effect of the electron-loop 
vacuum-polarization on a single Coulomb photon can be 
expressed as an integral over the vacuum-polarization 
function H(q 2 ). Here q may be either space-like or time- 
like. The first choice (q 2 < 0) leads to the integral 



A^E 



rf 3 



<l 



(27T) 



p(a 2 



[-n(-a] 



(3) 



Here p is equal to /5 2 p — P2S, f>2P and p 2 s being Fourier 
transforms of squares of non-relativistic Coulomb wave 
functions for the 2P and 25 states: 



P2P(2S) 



d S r\(j) 2P{2S) {r)\ 2 e 



iq-r 



Carrying out the integration we obtain 
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l-3a 2 



2a 4 



P2P 



(1 



,2W 



P2S = 



(1 + a 2 ) 



2\4 



(4) 



(5) 



where a = \q\/(Zam r ) and p 2 p is averaged over three 
degenerate states. 

The second choice (q 2 > 0) gives rise to the integral M 



/>OG 

A (II) E = m r (Za) 2 / dtu(t) 

J a 2(1 



2 



where 



and 



= 



m r a 



1 



0.737 383 76 (30) 



u(t) = -ImR{q A 



tm„ 



(6) 



(7) 



(8) 



Although Eqs. @ and @ are analytically equivalent, 
they are totally different from the viewpoint of numerical 
integration. Thus they provide a useful check whenever 
both real and imaginary parts of n are available. For di- 
agrams containing several vacuum-polarization loops in 
one Coulomb photon line, Eqs. (||) and (||) must be mod- 
ified accordingly. Insertion of vacuum polarization loops 



1 



in several Coulomb photon lines can be handled by the 
non-relativistic bound-state perturbation theory. 



FIG. 1. Three second-order vacuum-polarization dia- 
grams inserted in the Coulomb photon line exchanged by the 
muon and the proton. 

Let us hrst consider insertion of three second-order 
vacuum-polarizations in a Coulomb photon (see Fig.[l]). 
The contribution U (p2:3 \q 2 ) of this improper diagram 
can be expressed in terms of the second-order vacuum- 
polarization function Il( 2 )(g 2 ) as 



n^V) = (n< 2 V)) 3 



(9) 



where lT^ 2 ) is known analytically and has the spectral 
function 



. la 4m 2 / 2m 2 \ 9 „ 9 . 

uV(t) = 3-^1 + ? 2 >4™ 2 . (10) 

Substituting H^ 2 ^ in Eq. (g) and evaluating the integral 
numerically^] we find 

A (I) E (p2:3) = q 0Q6 253 4 ^ mr ( Za f _ (n) 

The result of the second method @ agrees with (|TT|): 
A (ii) E ( P 2:3) = q QQg 253 9 (10) m r (Za) 2 (^) 3 . (12) 

Another evaluation of A^E^ p2:3 ^ using the parametric- 
integral form of Il( 2 ) given in Ref. p] leads to 

A (/)£( P 2:3) = 0Q6 253 g (8) mr(Za) 2 ^ 3 _ (13) 





FIG. 2. Insertion of one second- and one fourth-order vac- 
uum-polarization diagrams in the Coulomb photon line ex- 
changed by the muon and the proton. 

The next contribution comes from diagrams involv- 
ing one second-order and one fourth-order vacuum- 
polarization insertions (see Fig. |^). This contribution 
is given in terms of 



n 



{pip2)( n 2\ _ 



(^) = -2n( 2 V)n (4 V) 



x This and subsequent integrals are evaluated numerically ei- 
ther on DECq or on Fujitsu- VX of NWU, or on both, by the 
adaptive-iterative Monte-Carlo subroutine VEGAS |H) . 



where lT^ 4 ) is the fourth-order vacuum-polarization func- 
tion §. Substituting Il( 4 ) into Eqs. (§) and © we obtain 

A (I) E ( P lp2) = 04 g 248 ( 5 ) mr ( Za )2 ^ 3 ^ (15 ^ 



and 



A (ii) E (p4 P 2) = Q46 243 ( 16 ) 



(16) 



We also evaluated A^E^p^ using the parametric- 
integral form of n^ 4 ) o: 



A (I) E (pip2) = QM6 25 q 



(2)m r (Za) 2 ^y . (17) 




FIG. 3. Sixth-order vacuum-polarization diagrams with a 
second-order vacuum-polarization inserted in the fourth-order 
vacuum-polarization diagrams. 

The third contribution comes from the sixth-order 
vacuum-polarization term JJ^p 4 ^ 2 )) obtained by inserting 
a second-order vacuum-polarization loop in the fourth- 
order vacuum-polarization diagram (see Fig.[|). The form 
of n( p4 ( p2 )) convenient for numerical integration is as in- 
tegral over Feynman parameters This can be done 
easily by adapting to the Lamb shift the program written 
previously for the electron g — 2 || . This leads to 

A (i) E ( P 4(p2)) = 013 628 ^ TOr ( Za )2 ( 18 ) 

The MS renormalized imaginary part of n( p4 ( p2 )) is 
known in a two dimensional integral form |{J . Converting 
it to the on-shell renormalized one and using Eq. (^) , we 
obtained 



A (n) E (pi(p2)) = 013 626 (1) mr ( Za f fSL 



(19) 





FIG. 4. Sixth-order vacuum-polarization diagrams with a 
single electron loop. 
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The fourth contribution comes from the sixth-order 
vacuum-polarization diagrams with a single electron 
loop. The exact form of this contribution is known only in 
a parametric-integral form Its imaginary part is not 
available in a form convenient for numerical work. We 
have therefore evaluated it using Eq. (|) only. There are 
eight topologically distinct diagrams (see Figji]). Each 
diagram can be written as a sum of various divergent 
terms and a finite part AH^ 6l \ where i = a,b, h. Af- 
ter renormalization the sum of these diagrams is free from 
any divergence and can be written as Q 

n( p6 > = 2(An (6a) + AIT (6c) + AIT (6d) + An (6/) ) 
+ An (6fc) + 4An (6e) + An (6s) + An (6/l) 

- 4AB 2 n (4) - 2[AB 4a + AL 4x + 2AL 4c 
+ AB ib + AUi + 2AL 4s + ^-(AB 2 ) 2 }U^ 



2(A6m 4a + A<5m 4h )n (2 * ) 



(20) 



where AB 2l are finite parts of renormalization 

constants and II^ 2 ) and II^ 4 ) are renormalized vacuum- 
polarization functions of second- and fourth-order, re- 
spectively. n^ 2 *) is the second-order vacuum-polarization 
function with a mass insertion vertex. Precise definitions 
of these functions are given in Ref. jl0| . The numerical 
values of the coefficients of I±( 4 ) , II^ 2 ) and I^ 2 *' are 



AS, = 



3 a 

4 7T 



AB 



^(AB 2 ) 2 = 0.871 680 (27) f- 

2 \ 7T 



4« 



A8m ia + ASm ib = 1.906 340 (21) 1-1 , (21) 

where the last two are new evaluations. The Lamb Shift 
contributions from HW, Il( 2 \ and n^ 2 *' can be easily 
obtained by numerical integration: 



AE {p4) = 0.045 922 7 (4) m r {Zaf 
AE {p2) = 0.017 452 8 (3) m r (Za) 2 



ae (p2*) = _ .009 001 8 (2) m r (Za) 2 -. (22 



The Lamb Shift contributions AE^ Qa \ coming 
from the ultraviolet- and infrared-finite parts of diagrams 
AIl( 6a ), • • •, are numerically evaluated. The results are 
summarized in Table I. The second and third columns 
list the results of integration carried out in double preci- 
sion and quadruple precision, respectively. The purpose 
of the latter calculation is to see whether the former indi- 
cates sign of losing significant digits due to rounding-off, 
which we call digit-deficiency problem and is the major 
source of uncertainty of on-the- computer renormalization 
J[jJ . The excellent agreement between two calculations 
shows that the estimated error of the former is not sig- 
nificantly affected by the digit-deficiency problem and 



can be safely assumed to be mostly statistical. We there- 
fore choose the double precision value, which has higher 
statistics, as our best estimate: 



A (i) E (pd) = o_oi7 410 (9) m r (Za) 2 (- 



(23) 



As a cross-check, we also evaluated AE^ p6 ^ using the 
Pade-approximation of the vacuum-polarization function 
from Ref. Q. We did this using both methods I and 
II. The [2/3] and [3/2] Pade approximations give nearly 
identical results. Taking their average we obtain 



AiI)E Pade = 0-017 414 9 (25) m r (Za) 2 (-] 
A(II)E Pade = 0- 017 414 9 ( 26 ) m r (Za) 2 (-Y . (24) 



These results are consistent with each other and agree 
with (p3| ) to three significant digits, or within one stan- 
dard deviation of (|22| ) . Obviously either ( p3| ) or ( pi[ ) has 
sufficient precision as far as comparison with experiment 
is concerned. Note, however, that the uncertainties given 
in are those resulting from numerical treatment of 
the Pade approximation and do not include those caused 
by the Pade method itself. It is argued in a separate 
paper |llj that the uncertainty of the Pade model itself 
is about 0.001 percent and hence the true value will be 
found well within the uncertainties given in (p4|). 




a 



FIG. 5. Representative sixth-order diagrams in which vac- 
uum-polarization insertion occurs in two and three Coulomb 
photon lines. Vertical lines represent the muon moving in the 
Coulomb potential generated by the proton, which is indi- 
cated by a "x". 

Thus far we considered only diagrams in which one 
Coulomb photon line is modified by the electron-loop vac- 
uum polarization. Additional contributions of order a 3 
arise from the diagrams of Fig. ^| in which two and three 
Coulomb photons are modified by vacuum polarization. 
Their contributions to the Lamb shift can be found by 
the bound-state perturbation theory: 

A£(Fig.5a) = 0.009 166 (2) m r (Za) 2 (-) , 
A£(Fig.5b) = 0.024 805 (3) m r (Za) 2 , 
A£(Fig.5c) = 0.002 535 (1) m r (Za) 2 ■ (25) 
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In this calculation we used the reduced non-relativistic 
Coulomb Green function for 2S and 2P states given by 
Eqs. (23) and (24) of Ref. §. 

Collecting @, @, (||), (||), and (fl), we obtain 
the total contribution to the Lamb shift (R ) due to the 
sixth-order vacuum-polarization effect. 

Evaluation of various lower-order contributions to the 
2-Pi/2 — 2S*i/2 Lamb shift C of the muonic hydrogen are 
summarized in Ref. J^Q. In addition we have obtained 
the hadronic vacuum-polarization correction of 0.0113(3) 
meV following Ref. JT2J . These results and our result (|l|) 
lead to the most precise theoretical prediction 



C = (206.068 (2) - 5.197 5 rl) meV, 



(26) 



where r p is the proton charge radius in units of fm. The 
uncertainty in the first term of ( ^6| ) is our estimate of 
theoretical error. 

To improve the theoretical prediction further, it is nec- 
essary to have better estimate of the effect to the Lamb 
shift and hyperfine structure of the muonic hydrogen due 
to the proton's internal structure beyond elastic form fac- 
tors. Recently the proton polarizability correction to the 
hyperfine structure of the hydrogen and muonic hydrogen 
was obtained |l3|. There are also references for ordinary 
hydrogen and deuterium [ ^i|JT^ ] . Unfortunately they are 
not directly applicable to the muonic hydrogen because 
of very different energy scale. 

Measurement of C to 10 ppm, or 0.002 meV, will lead 
to improvement in the value of by an order of mag- 
nitude over those determined from the elastic scattering 
form factor measurements, making it possible to resolve 
the long-standing discrepancy between |l6| and |l7j . The 
new value of will also play an important role in testing 
the validity of QED in terms of high precision measure- 
ments of the hydrogen atom jl8|]. Another impact of 
accurate determination of will be to stimulate evalua- 
tion of Tp from the lattice QCD more precise and reliable 
than those available at present |l9| . 
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TABLE I. Contributions to the 2P 1 / 2 — 2Si/2 muonic hy- 
drogen Lamb shift from the sixth-order vacuum polarization 
diagrams with a single electron loop. The overall factor 
772. r (Zafia/nf is omitted. The second and third columns 
give results of integration in double precision and quadruple 
precision, respectively. Their difference is listed in column 4. 



Term 


Doub. precis. 


Quad, precis. 


Difference 


AE {6a) 


0.044 769 (4) 


0.044 739 (51) 


0.000 030 (52) 


AE {6b) 


0.028 654 (4) 


0.028 640 (35) 


0.000 014 (36) 


AE {6c> 


-0.025 393 (3) 


-0.025 368 (23) 


- 0.000 025 (24) 


AE {6d> 


-0.026 376 (2) 


-0.026 371 (21) 


- 0.000 005 (22) 


AE ((ie> 


0.151 356 (4) 


0.151 334 (46) 


0.000 022 (47) 


AEW> 


-0.067 139 (3) 


-0.067 144 (30) 


0.000 005 (31) 




0.019 536 (3) 


0.019 540 (23) 


- 0.000 004 (24) 


AE (6h> 


0.025 877 (2) 


0.025 858 (22) 


0.000 019 (23) 



2 K. Pachucki informed us that F. Kottman pointed out that 
the sum of all contributions listed in Ref. Q was 206.049 meV, 
not 205.932 meV. 
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